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ABSTRACT 

A general formalism is developed for calculating the luminosity function and the expected number of ob- 
served GRBs above a peak photon flux S for an arbitrary GRB jet structure. This formalism directly provides 
the true GRB rate for any jet model, instead of first calculating the GRB rate assuming isotropic emission and 
then introducing a 'correction factor' to account for effects of the GRB jet structure, as was done in previous 
works. We apply it to the uniform jet (UJ) and universal structured jet (US J) models for the structure of GRB 
jets and perform fits to the observed log A^- log 5 distribution from the GUSBAD catalog which contains 2204 
BATSE bursts. We allow for a scatter in the peak luminosity L for a given jet half-opening angle 6j (viewing 
angle 6'obs) in the UJ (USJ) model, which is implied by observations. A core angle 9c and an outer edge at 6'max 
are introduced for the structured jet, and a finite range of opening angles 9^^^ < 9j < 6',nax is assumed for the 
uniform jets. The efficiency for producing 7-rays, e~^, and the energy per solid angle in the jet, e, are allowed 
to vary with 9j (6'obs) in the UJ (USJ) model, e-y cx 6*"* and e cx 9~". We find that a single power-law luminosity 
function provides a good fit to the data. Such a luminosity function arises naturally in the USJ model, while 
in the UJ model it implies a power-law probability distribution for 9/, Pifij) cx The value of q cannot be 
directly determined from the fit to the observed log A^- log 5 distribution, and an additional assumption on the 
value of a or b is required. Alternatively, an independent estimate of the true GRB rate would enable one to 
determine a, b and q. The implied values of 9c (or ^min) and 0^^^, are close to the current observational limits. 
The true GRB rate for the USJ model is found to be «grb(z = 0) = 0.86![| Gpc"^ yr"' (1 a) while for the UJ 
model it is higher by a factor f{q) which strongly depends on the unknown value of q. 

Subject headings: gamma rays: bursts — gamma-rays: theory — cosmology: observations — ISM: jets and 
outflows 



1. INTRODUCTION 

There are several lines of evidence in favor of jets in 
gamma-ray bursts (GRBs). For GRBs with known redshift 
z, the fluence can be used to determine the total energy output 
in 7-rays assuming spherical symmetry, /s^jso- The values of 
ii-yjso that were inferred in this way sometimes approached, 
and in one case (GRB 991023) even exceeded Mqc^. Such 
high energies are very hard to reconcile with progenitor mod- 
els involving stellar mass objects. A non-spherical, colli- 
mated outflow (i.e. a jet) can significantly reduce the total 
energy output in 7-rays compared to /s-y.iso, since in this case 
the 7-rays are emitted only into a small fraction, //, ^ 1, of 
the total solid angle. A more direct (and probably the best 
so far) line of evidence in favor of jets in GRBs is from 
achromatic breaks in the afterglow light curves jRhoadsl 19971 
[1999; Sari. Piran & Halpern 1999). 

Despite the large progress in GRB research since the dis- 
covery of the afterglow emission in early 1997, the structure 
of GRB jets is still an open question. This is a particularly 
interesting and important question, since it affects the total 
energy output and event rate of GRBs, as well as the require- 
ments from the central engine that accelerates and colhmates 
these relativistic jets. 

The leading models for th e jet structu re are: (1) 
the uniform jet (UJ) model JGranot et alJ 1200 IL l2002t 
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Kumar & Panaitescul 120001: iModerski, Sikora & Bulik 200C ; 
Panaitescu & Mf eszfarodl J[T999i: .Rhoads, .1997. .199S : 
Sari. Piran & Halpern 1993), where the energy per solid 
angle e and the initial Lorentz factor Fq are uniform within 
some finite half-opening angle Oj, and sharply drop outside 
of 9i, and (2) the universal structured jet (USJ) model 
(^Lipunov. Postnov & Prokhorov'200l HRossi. Lazzati & ReesI 
2002; Zhana & Meszaros 2002), where e and Fq vary 
smoothly with the angle 9 from the jet symmetry axis. In the 
UJ model the different values of the jet break time tj in the 
afterglow light curve arise mainly due to different 6j (and to 
a lesser extent due to different ambient densities). In the USJ 
model, all GRB jets are intrinsically identical, and the differ- 
ent values of tj arise mainly due to different viewing angles 
6'obs from the jet axis [in fact, the expression for tj is similar 
to that for a uniform jet with e ej Oobs) and 9j 6'obs]- The 
observed correlation f, oc oom. Frail & Kulkarnil 

2003; Frail et al. 2001) implies a roughly constant true 
energy E between different GRB jets in the UJ model, and 
e (X 9~^ outside of some core angle 9c in the USJ model 
jRossi. Lazzati & Rees 2002; Zhang & Meszaros 2002).^ 
The probability distribution of jet half-opening angles be- 
tween different GRBs in the UJ model, P{9j), is not given 
a-priori by the model, and is a free function to be determined 
by observations. 

Several methods have been used so far in order to constrain 
the structure of GRB jets and help distinguish between the UJ 
model and the USJ model. The afterglow light curves are 
similar for the UJ and USJ models, but nevertheless some 

^ The latter is obtained assuming that the efficiency in producing 7-rays, 
e-y, does not depend on 8. We later examine the consequences of relaxing 
this assumption. 
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differences still exist which m ight help distinguish between 
them iGranot & Kumail2003t fKumar & Granot 2003). After- 
glow light curves also constrain the jet structure in the USJ 
model, as discussed in |3l This method requires a very good 
and dense monitoring of the afterglow light curves, especially 
near the jet break time tj. Another possible way to distinguish 
between these two models for the jet structure is through the 
different expected evolution of the linear polarization during 
the afterglow (Rossi et al. 2004). However, some difficulties 
and complications exist in this method, like the poor quality 
of most polarization light curves, and a possible ordered mag- 
netic field component that might cause polarization that is not 
related to the jet structure JGranot & Koni^ll|2003 ). The dis- 
tribution of viewing angles 9ohs, inferred from the observed 
values of f ,, has been used to argue in favor of the USJ model 
(|Perna, Sari & Frail 2003). However, when the known red- 
shifts of the same sample of GRBs are also taken into ac- 
count, there is a very poor agreement with the predictions of 
the USJ model (Nakar. Granot & Guetta 2004). In order to 
reach strong conclusions using this method, a large and uni- 
form sample of GRBs with known redshifts is needed, like the 
one expected from Swift. 

In this paper we use the observed log A^- log 5 distribution 
(the number of GRBs above a limiting peak photon flux 
S) in order to constrain the jet structure. A somewhat sim- 
ilar analysis was done by Guetta. Piran & Waxman (20041. 
They found that the observed log A^- log 5 distribution rules 
out the USJ model, due to the paucity of GRBs with small 
peak fluxes compared to the prediction of the USJ model, 
and fit a double power-law luminosity function for the UJ 
model. Firmani et al. (2004) tried to constrain the redshift 
evolution of the GRB rate and luminosity function, using the 
logA^-log5 distribution (as well as the redshift distribution 
derived from the luminosity-variability relation). 

This paper improves upon previous works by: (1) allow- 
ing the efficiency in producing 7-rays to be a function of the 
angle 6*, Cy = e^( 9), and varying e{9) accordingly, using the 
iFrail et aP ( 120011) relation, (2) including an internal dispersion 
in the peak isotropic equivalent luminosity L at any given an- 
gle 9, (3) introducing an inner core angle 9^ and an outer edge 
9m!ix in the USJ model, (4) using a larger and more uniform 
GRB sample. 

In ^we allow e-y and e to vary with 9, and derive the 
constraints on the power-law indexes that are implied by the 
IFrail etalJ(l2()0r ) relation. The formalism for calculating the 
luminosity function and the observed GRB rate as a function 
of the limiting flux for different jet structures is derived in ^ 
In l^and ^we compare the observed logA^-log5 distribu- 
tion to the predictions of the USJ and UJ models, respectively. 
Our results are discussed in ^ 

2. THE ENERGY AND GAMMA-RAY EFFICIENCY DISTRIBUTIONS 

The efficiency for producing 7-rays, e-y, is taken to be a 
function of the angle 9 from the jet symmetry axis in the 
USJ model, and a function of the jet half-opening angle 9j 
in the UJ model. For convenience we assume a power- 
law dependence on 9 in the USJ model, £-^(9) = Q{9^^^ - 
9)ey omm[l ,{9 /9c)~''^^'], where Q(x) is the Heaviside step 
function. An outer edge at ^^ax has been introduced, as weU 
as a core angle, ^ 1 /rmax ~ 5 x 10""* rad, which is needed 
in order to avoid a divergence at = 0. Here r^ax ~ 2000 is 
the maximum value of t he Lorentz factor to which the fire- 
ball can be accelerated (.Guetta. Soada & Waxmaa2001) . For 
the UJ model P(9j) is restricted to a finite range of values. 



^^min < < ^max, and in analogy to the USJ model we chose 

eyi9j) = e(^„,ax - 9j)e{9j - 9^i„)eyfl{9j/9^,,r''"' ■ 

The energy per solid angle is also assumed to behave as a 
power law, e(9) = Q{9^^^-9)eormn[l,{9/9cT"''^'] in the USJ 
model and e(9) = e(9^,,-9j)e(9j-9,^,,)eo (9j/9^i„r'"' in 
the UJ model. The power-law indexes a and b can be differ- 
ent between the USJ and UJ models (as is emphasized by the 
different subscript for the two models), and we also consider 
different values for a and b within each model. 

The external mass density is taken to be a power-law in the 
distance R from the central source, pext = AR~''. For a con- 
stant efficiency jb = 0) the observed correlatio n tj cx E~\^^ 
('Bloom, Frail & Kulkarni 2003; Frail et aVlOOl) implies a = 
2 (Rossi, Lazzati & Rees 2002; Zhang & Meszaros 2002). In 
the foll owing, we allow b to vary and find the joint constraint 
that the lFrail et alJ ( 1200 Ih relation puts on a and b. 

The afterglow emission becomes prominent when the GRB 
ejecta sweeps enough of the external medium to decelerate 
significantly. After this time most of the energy is in the 
shocked external medium, and energy conservation implies 



(1) 



where fi is the swept-up rest mass per soUd a ngle. 

In the UJ mod el T( ti) ~ l/9j dRhoadsl (ToOTI: 
|Sari, Piran & Halperd Il999h while in the USJ model 

T(tj) ^ l/9ob, (Rossi, Lazzati & Rees 2002; Zhang & 
Meszaros 2002) and the emission around tj is dominated by 
material close to our line of sight. Together with the relation 
we have 



1 

4^ 



{3-k)e(9) 



Q2{A-k)IO-k) ^ Qm-k)-a]ICi-k) ^ (2) 



where 9 = 9j for the UJ model and 9 = 6'obs for the USJ 
model. Now, in order to satisfy the observed correlation^, 
'""^'^ we must have 



X = a + b-- 



2i4-k)-a 
3^ 



(3) 



Therefore if we allow the efficiency and the energy per solid 
angle to vary in this way, we obtain the condition b = (2- 
a){4-k)/(3-k)oT a = 2-b(3-k)/{4- k) in o rder to reproduce 
the observed correlation of Frail et al. (2001). 

For the USJ model, interesting constraints on the power- 
law index ausj of the energy per solid angle h ave been derived 
from the shape of the afterglow light curve ( iGranot & Kumad 
[2903 ; Kumar & Granot 2003). For ousi ^1-5 the change in 
the temporal decay index across the jet break is too small com- 
pared to observations, while for ausi 2.5 there is a very 
pronounced flattening of the light curve before the jet break, 
which is not observed. The latter feature arises since the inner 
parts of the jet near the core, where e is the largest, become 
visible. This explains why this feature becomes more pro- 
nounced for flusj > 2, where most of the energy in the jet is 
concentrated at small angles, near the core. There is no sharp 
borderline where the light curves change abruptly. Instead, 
the light curves change smoothly with the parameter ausj- Al- 
together, the observed shape of the afterglow light curve con- 
strains the parameter ous: to be in the range 1 .5 < ausj ^2.5. 

* We have E-\^^ oc since E^^^^aid) = 4Trt{9)t^{e). The peak isotropic 
equivalent luminosity is given by L = E^ \^„/T where T is the effective dura- 
tion of the GRB, which is assumed to be uncorrelated with L so that L has the 
same 6-dependence as £-y.iso 
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It is important to stress that this constraint applies only to the 
USJ model and not to the UJ model, since in the former ausj 
determines the structure of an individual jet (and therefore af- 
fects the light curves) while in the latter the structure of an 
individual jet is fixed and ouj only affects how the uniform 
e changes between different jets of different half-opening an- 
gles Oj. 

One might also try to constrain the power-law index b of 
the gamma-ray efficiency from observations. This would, 
of course, apply to both the USJ model and the UJ model. 
The way in which this has been done so far is by taking 
e-y = £'-y4so/(£7,iso + £k,iso) whcre isk.iso is the initial value of 
the isotropic equivalent kinetic energy. However, it is hard 
to evaluate /Sk.iso very accurately. It is usually evaluated from 
afterglow observations several hours to days after the GRB, 
and can provide an estimate accurate to within a factor of ~ 2 
or so for the kinetic energy at that time. An additional and 
potentially larger uncertainty arises since at early times there 
is fast cooling and radiative losses might reduce the initial ki- 
netic energy by up to an order of magnitude or so. This can be 
taken into account, but introduces an additional uncertainty in 
th e value of isk isn that is estimated in this way. 

iPanaitescu & KurnaJ (1200 L .200 2) evaluated /ik.iso from a 
fit to the broadband afterglow data for ten different GRBs and 
obtained > 0.5 for most of these GRBs and e-y > 0.1 for 
all of them. There seems to be no particular correlation with 
9, but due to the small number of bursts and the reasonably 
large uncertainty in /Sk.iso there might still be some intrinsic 
correlation. Llovd-Ronning & Zhang (2004) used the X-ray 
luminosity at 10 hr to estimate the kinetic energy at that time 
and used a simple analytic expression (Sari 1997) to account 
for the radiative losses. In this way they estimated for 17 
GRBs and one X-ray flash (XRF 020903). Their Figure 6 
suggests that decreases with 9, although they claim that 
this is not statistically significant. Since most of the estimates 
for e-y are > 0.5 and are for small values of 9, and since e-y < 1 
by definition, this suggests b >Q (otherwise we would have 
e-y > 1 at large 9, which is impossible). Moreover, Figure 6 
of Llovd-Ronning & Zhang (2004) suggests that < < 1. If 
we assume that <b < 1, then Eq. (|3j would imply 1.25 < 
fl < 2, 2 < A < 2.25 for ;t = 0, and 1.5 < a < 2, 2 < A < 2.5 
for k = 2. 

3. LUMINOSITY FUNCTION FOR DIFFERENT JET STRUCTURES 

For simplicity, we assume that the emission from the GRB 
jet is axially symmetric, so that the observed luminosity de- 
pends only on the viewing angle 9 from the jet axis, and does 
not depend on the azimuthal angle ip. It is convenient to de- 
fine the probability distribution P{L\9), where P{L\9)dL is the 
probability for an isotropic equivalent luminosity between L 
and L + dL, when viewing the jet from an angle 9. Naturally, 
we must have J P(L\9)dL = 1 for any value of 9. The proba- 
bility of viewing the jet from an angle between 9 and 9 + d9 is 

simply P{9)d9 = sm9d9, and ^^'^P(9)d9= 1. Avera ging over 
the viewing angle 9 we obtain the overall probability distribu- 
tion for L, 



P(L): 



P(9)P(L\9)d9 = P(L\9) d cos 9 . (4) 



For a universal GRB jet structure, i.e., if all GRB jets have 
the same intrinsic properties, then Eq. (0} represents the GRB 
luminosity function. 
As a simple and very useful example, let us consider a 



universal structured jet model where^ ^"(^1^) = S[L-L{9)] 
and L(9) = Lo(9/9oT^. Using Eq. (gj this implies P(L) = 
(6lo/ALo)(L/Lo)"'"'/^ sin6'(L), where 9{L) = 9uiL/LoT^/^, and 
for 9(L) < 1 it reduces to P(L) « (6'g/ALo)(L/Lo)"'"^/^- Al- 
ternatively, one may assume^ P{L) = CqL'"^, where 77 = 1 +2/ A 
and Co ^ X~^9f^LQ~ is determined by the normalization con- 
dition, / P(L)dL= 1. 

An important point to stress here is that P(L) represents the 
average over all possible viewing angles. Even if the jet is 
assumed to have a sharp outer edge at some finite angle ^max, 
with some probabiUty distribution P*(L) for 9 < ^max, such 

that P*(L) = (1 -cos6lmax)"' f^"" d9sin9P(L\9) and PiL\9 > 
9j) = S(L), then P(L) = (l-cos9^^^)P*(L) + cos9^,J(L). 

If the GRB jet structure is not universal, and P{L\9) de- 
pends on additional parameters that describe the jet struc- 
ture, then Eq. ^ needs to be averaged over these parame- 
ters in order to obtain the GRB luminosity function. For ex- 
ample, for a uniform jet of half-opening angle 9j, we have 
P(L\9,9j) = Q(9j-9)P*(L\9j)+e(9-9j)S(L) where P*(L\9j) 
is the probability distribution of L inside the jet (at 9 < 9j) for 
a given value of 9j. In this case, if Pi9j) is the probability 
distribution for 9j, then the GRB luminosity function is given 
by 

P(L)= / P(9j)d9j / P(9)d9P(L\9,9j) 
Jo Jo 

P{9j)d9j [(l-cos6lj)P*(L|6i;) + cos6'/(L)] (5) 



In analogy with the USJ model, we assume for the UJ model 
that e = e{)(6lj/6'mi„)-"™ and = t^.o{9j/9^,^y'"'' . Note that 
since e = £'iso/47r « £'/27r0j this implies that the true energy 
is not necessarily constant, E oc 0^"°"' . In order to imitate a 
structured jet with a uniform core, we chose 
P{9j) = 0(9 J- 0™„)e(0„,x - 0j)C(q)9-/ + B(q)S(9j - 9^,) , 

(6) 

where C(q) = [l-B(?)](l-?)/(0,n7x-^i7D for ? ^ 1 and C(l) = 
[1 -B(l)]/ln(6'max/6'min) from the normaHzation / P(9j)d9j = 
1. For P*{L\9j) = S[L-L{9j)] where L{9j) = L^ini9j/9^,,r^"' 
and Lmin = i(6'max), this implies 



P(L) = e(L-L™„)e(L„,x-i) 



AlIlLmin 



[l-cos9j(L)] 



-l-(l-?)/Auj 



+ B(1 -COs6'n,in)(5(L-Lniax) 



Bcosl 



d9i 



S(L) 



ie(L-L^.n)e(L„ax-^) ^, — 



-\-0-q)IXm 



+B^5(L-L^,,) 



2 r (9^~'' -9^''') 



j_^_mm 

2 2 (3-^) 



6iL), 



(7) 



^ This is the standard universa l structured jet model IRossi. Lazzati & Reed 
I2002:lzhang & Meszaro J2002I) for A = 2, and assuming no scatter in L for a 
given 9. 

^ This is a pure power law in L for all 0, instead of just for S <C 1 as we 
had before. 
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P*{L) = Q{L-L^^)Q{L^^-L) 
+BS(L-L^^^) , 



Auji'min 



-l-(l-?)/Auj 



(8) 

where 9j(L) = 0„ax(i/im,„)"'/^™. 

In the case of the structured jet discussed above, with a 
sharp outer edge at some finite angle 0max, and a uniform core 
within some angle dc, we have 

P(L\e) = Q(0,- 0)S(L - L™, J + 8(0 - e^,,)6(L) 

+Q(9-0,)Q{9^,,-9)6[L-Li0)], (9) 

where L(9) = Lmin(6'/6'max)"^"'' and L^in = L(0max)- Therefore, 



P{L) = e(L, 
+ (1 

~ G)(^max 



IX 

COs6'c)(5(L-Lniax 

L)e(L-L,„,„) 



fi'max sin 6'(L) 



AuSjimin 

) + (5(L)cos6'n 



-I-l/Ausi 



92 

^max 



.)+ 1 



AuSJ^rnin 
"max 



-1-2/Au 



(10) 



P*{L): 



Q{L^^-L)Q(L-L^,) 



(1- 



-1 

'max' 



sin 6'(L) 



-1-1/Au 



(l-cos6'c) 



(l-COS^max) 

'.Q{L^,^-L)Q{L-U 



^{L i'max) 



IT 



-1-2/Au 



5(Z^ Z^niax) ; 



(11) 



where 9(L) = 6'max(i/imin)"'/^"" and P*(L)dL = 1 as it 
should be. Note that the coefficient in the first term in the ex- 
pression for P{L) includes fl^ax only through the combination 
^'max^i'';^ = 6'max[i'(6'max)] which is independent of 9^^^ (for 
a fixed normalization of L) since L{9) cx 6*"^. 

By comparing equations Q and (I10> one can also see that a 
similar luminosity function can be obtained for different val- 
ues of A = « + as long as Auj/Ausj = (3-q')/2. In this case, 
when q^l then Auj 7^ Ausj- For the same ^^ax, L^in and L^ax 
we have Ln,ax/imi„ = (6'max/6'mi„)^™ = (6'max/6'c)^"'' which im- 
pUes 9c/9^in = (6'min/6'max)*'"'''''^ (whcrc 6*, < 6ln,i„ for ^ < 1 
and dc > 6'n,in for q> 1). In order for the luminosity functions 
to be the same also for the core of the structured jet, which is 
represented by the term cx (5(L-Lniax), the ratio of this term to 
the other terms should be the same for equations Q and (I10> . 
Thus we obtain 



C(q)-- 



min 



qI-9 _ /5I-9 
7max 'J„iin 



B(q)=<l + 



0-q) (l-q) 



(l-q) 



(12) 

(13) 



whereS(l) = C(l)/2= [l+21n(0„ax/^min)]"' and C(^)/B(^) = 

(3-?)Cr „ 
The normalization of the first two terms in Eq. (Q for the 

UJ model is smaller than that of the (corresponding) first two 



terms in Eq. ilQ\ for the USJ model by a factor of 



9 

-'max 



l-q 



0-q) 
(l-q) 
(3-q)/(l-q) 



-'mill 

9 

'max 



l-q- 



(q<l), 



.[2/(^-l)](^,„ax/^mi„)''-' (^> 1) 



(14) 



where /(I) = 1 +21n(6'inax/^min)- This factor / is the ratio of 
the true GRB rate for the UJ model to the true GRB rate for the 
USJ model that corresponds to the same luminosity function 
(i.e. Eos. I71andll0>. Figure ^ shows /(^) for three different 
values of 6',nax/ 9 min- The fact that f(q) > 1 means that a larger 
number of uniform jets is needed, compared to structured jets, 
in order to reproduce the same luminosity function and the 
same observed rate of GRBs, Ncrb- This means that for the 
same A'grb the intrinsic GRB rate per unit comoving volume, 
Rgrb(z), and specifically Rgrb(z = 0), is larger by a factor of 
/ for the UJ model compared to the USJ model. 

The same factor / is obtained from the ratio of the energy 
output in gamma-rays of a single structured jet to the average 
energy output in gamma-rays of a uniform jet. This ratio must 
equal the ratio of the intrinsic rates since the total energy out- 
put in gamm a-rays pe r unit time per u nit volume must be the 
same dGuett a. Piran & Waxmanl2004 . 

We can see that / '--^ 1 for < 1 and / ^ (^max/^'min)*"' 
for q > I. This can be understood as follows. Since most 
of the solid angle of the structured jet is near 9^^%, the num- 
ber of uniform jets with 9j ^ 6'max should be comparable to 
the number of structured jets. Therefore, for q < I where 
most of the uniform jets have 9j 9max we have / 1, while 
for q > I where most of the uniform jets have 9j ^ 9mm, 
f ^ (fi'raax/^min)'' ' IS roughly the inverse of the fraction of 
jets with 9j ^ 6'max- 

We have demonstrated that the luminosity function of the 
USJ model can be imitated by the UJ model with the appro- 
priate choice of P(9j). The implied true GRB rate Rgrs(z) 
for the same observed GRB rate Nqrb would always be larger 
(by a factor /) for the UJ model. For the UJ model, the term 
B5{9j-9mm) in P(9j) which produces the term cx d(L-Lmax) in 
P{L) is somewhat artificial, and was introduced only to make 




Fig. 1 .— The ratio / (Eq. [141 of the true GRB rate for the UJ model to 
the true GRB rate for the USJ model, given the same observed GRB rate and 
luminosity function (i.e., the same logA? — logS distribution), as a function 
of the slope q of the probability distribution for the opening angle dj of the 
uniform jet, P(d j) oc d^'' for Sniin < Sj < Bmux- See text for details. 
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a complete analogy with P(L) for the US J model. The exact 
form of the cutoff near 6'min is not clear and probably would 
not have a very large effect on the logA^-log5 distribution, 
making it hard to distinguish between the two models in this 
way. It is also important to keep in mind that the observational 
constraints on the parameter a from the shape of the afterglow 
light curve, namely 1.5 < ausj ^ 2.5, apply only to the USJ 
model and not to the UJ model. 

The form of the luminosity function given in equations (7) 
or ilQ\ is valid only if 9^L{6) ^ const, where 6 = Oobs for the 
USJ model and 9 = Oj for the UJ model. How well this con- 
dition is satisfied may be tested by using the values of both 
L and 6 that were estimated from observations for a sam- 
ple ofGRBs^_\^_considera sub-sample of 19 GRBs out of 
the IB loom. Frail & K ulkarm (2003) sample, for which there 
is both a known redshift and an estimate for 6. There are 
larger uncertainties in the determination of the peak lumi- 
nosity,^ as the bursts were detected by different instruments 
wit h different temporal and spectr al sensitivities. Follow- 
ing Guetta. Piran & Waxman ( 2004) we have extrapolated the 
ORB fluxes to the BATSE range (50keV-300keV) using the 
method described in Sethi & Bharaavi ( 2001 ). In our analysis 
we have used the median value of the spectral photon index in 
the 50 keV - 300 keV band for the long bursts sample, -1.6, 
as found by 'Schinid? 1*200 1). The redshifts and fluxes were 
taken from the table given in Van Putten & Reaimbau (2003i). 

We calculate the values of O^L for such a sample and for 
different values of A within the range allowed for the USJ 
model from the constraints given above (namely, 1 .5 < ausj ^ 
2.5 and Eq. |3}: 1 .83 < A < 2. 17 for A: = 0, and 1 .5 < A < 2.5 
for k = 2. Therefore overall 1.5 < A < 2.5. The result of this 
analysis is that the distribution of 9^L(9) is not quite a delta 
function, and there is some dispersion around the mean value. 
This dispersion is reasonably fit by a log-normal distribution. 



P(L\9): 



1 



Lav In 



X exp ' 



[lnL-lnL(6')]^ 
2^2 



(15) 



with L{9) = Lq{9/9q) ^, where 9qL() is the average value of 
9^L{9) for the observed sample, and a is the standard devia- 
tion of \vi{9^L{9)], which is determined by a fit to the disper- 
sion of the observed sample. This distribution approaches a 
delta function (i.e. with no scatter) for a ^Q. 

We perfo rmed fits of 9'^L to a sample of 19 bursts from 
iBloom. F^ il & Kulkarni (2003) which were observed by 
BATSE and Beppo-SAX. For A = (1.5,2,2.5,3) we ob- 
tained 9^^U_) = (15,5.8,2.2, 1.1) X lO^** erg s"' and a^bs. = 
(1.2, 1.2, 1.3, 1.5). These values of CTobs can be thought of as 
upper limits on a (i.e. a < Cobs), since the observed scatter in 
9^L{9) includes both the intrinsic scatter (which should be re- 
flected in a) and additional scatter due to measurement errors 
in 9 (which produces some scatter (Ten), which can be as high 
as tens of percent. Blo om. Frail & Kulkarni (2003) obtain a 
factor of 2.2 for the scatter in efl^, implying that (j\ne ^ 0.4, 
and therefore tJeiT S 0.4A. Assuming that the measure- 
ment errors in 9 are uncorrelated with the intrinsic scatter 



in 9 L(9), this implies 



2 

obs 



o- + CTe„ < + (0.4A)2and 



therefore (1.0, 0.9, 0.8, 0.9)<CT<(1. 2, 1.2, 1.3, 1.5) for A = 
(1.5,2,2.5,3). 

' We use the peak luminosity since it is usually the most relevant quantity 
for the triggering of the different detectors. 

This is since the dominant measurement error in comes from 9^, 
while L is measured more accurately. 



The luminosity function in this case can be obtained using 
equations ^ and (I15> . It is important to note that the integral 
in Eq. (|4} must be done over all possible viewing angles, even 
if the jet is assumed to have a core angle, 9^, and maximal an- 
gle, (?max- In fact, together with the power-law index (A) and 
normalization {9qLo), the other two intrinsic parameters that 
define a power-law universal structured jet are the angles of its 
outer edge (0max) and inner core {9c). For the UJ model out- 
lined above, 9^^^^ and 9^ are replaced by 6'tnax and 6'min which 
determine the range of possible 9j values, where P(9i) oc 9^ 
also depends on an additional parameter q. 

The observed rate of GRBs (over the entire sky) with a peak 
photon flux greater than S is given by: 



A^grb(>5)= / P(L)dL 



Rcrb(z) dV(z) 



dz , (16) 



70 (l+z) dz 

where P(L) is given by equations ^ and il5\ . z is the redshift 
of the GRB, Zmax(L,S) is the maximal redshift from which a 
GRB with luminosity L and peak flux S can be detected, and 
Rcrb{z) is the (true) GRB rate per unit comoving volume V . " 
The factor ( 1 +z)~' accounts for the cosmological time dilation 
and dV(z)/dz is the comoving volume element. 

4. CONSTRAINTS ON THE LUMINOSITY FUNCTION OF THE USI 

MODEL 

We consider all the GRBs in the GUSBAD catalog 
('Schniidt; 2004), which lists 2204 GRBs detected at a time 
scale of 1024 ms. This catalog contains all the long GRBs 
(Tgo > 2 sec; [Kouveliotou et al. 1993), detected while the 
BATSE on-board trigger ( Paciesas et al. 1999) was set for 5.5 
(7 over background in at least two detectors, in the energy 
range 50 - 300 keV. Using this sample we estimate the ra- 
tio Cmax/Cmin for cach burst, where Cmax is the count rate in 
the second brightest illuminated detector and Cmin is the min- 
imum detectable rate. We find (V/Vmax) = 0.335 ± 0.007. 

We consider also th e Rowan-Robinson SFR (RR-SFR; 
iRowan-RobinsonI 1 9991) that can be fitted with the expression 

Rcrb(z) = PO X ■{ ^^0 75 (17) 



z> 1 



where po = Rcrb(z = 0) is the true GRB rate per unit comoving 
volume at z=0. Throughout the paper we use the following 
cosmological parameters: (nM,^A,h) = (0.3,0.7,0.7). 

Objects with luminosity L observed by BATSE with a flux 
limit Slim are detectable to a maximum redshift Zmax(i','S'iim)- 
The limiting flux has a distribution G(5ii,n) that can be ob- 
tained from the distribution of C,ni„ of the GUSBAD cata- 
log. Considering five main representative intervals we ob- 
tain that 30%, 30%, 10%, 20%, and 10% of the sample have 
5iim - 0.23, 0.25, 0.22, 0.26, and 0.27 phcm'^ s"', respec- 
tively. Therefore, we have 



P{L)dL / G{Siim)dSiim 



A^grbO S) = 







P(L)dL 



JO 



(l+z) dz 
''"^ Rgrb(z) dV(z) 



dz 



(l+z) dz 



G{Slim)dSlim 



max(Z.,Sli„,) 



RcRBiz) dV(z) 
(l + z) dz 



dz 



dz 



,(18) 



' ' Note that Rqrb (z) is the true GRB rate, and it is obtained directly from 
this formalism without any need for a 'connection factor' etc. 



6 



Table 1. Fits to the logW-logS Distribution 



P{L) 
model 


model 
parameters 


Rgrb(z = 0) 
[Gpc-3 yr-1] 


goodness 
of fit 


single 
power 
law 


Ausj = 2.9^l,\ 
a = 0.5:°;? 
0,=O.O5+oof 
emax=0.7±0.2 


0.86:0;'^ (USJ) 
5.4+|8(UJ,^=l) 

— l.U ' J ' 


= 11.05 
(9 d.o.f) 
/" = 0.723 
(1.10 0-) 


single 
P.L.+ 
fixed 
Ausj 


Ausj = 2 (fixed) 
a = 0.8+{;'^ 
e,=0.03+00f 
0max= 0.5 ±0.05 


0.75:0;0^ (USJ) 
5.0:i;3 (UJ, q = 1) 


= 14.67 
(10d.o.f) 
P = 0.855 
(1.46 a) 


double 
power 
law + 
cr = 


o = 0.6±0.1 
/3 = 0.8^-2 


10.7 ±2.3 
(UJ, q=l) 


= 10.57 
(10d.o.f) 
P = 0.608 
(0.86 a) 



Note. — The best fit parameters with their 1 a confidence intervals are 
shown together with the implied true GRB rate and the goodness of fit. The 
confidence intervals are the projection onto the relevant parameter axis of the 
region in the multi-dimensional parameter space around the global minimum 
X?„i„ of where Ax^ = X^-xl^in < 1- 

tHereLji =L,/(105l ergs"'). 



where P{L) is given by equation (0}. For the USJ model 
P{L\e > 6l„ax) = S(L) while P{L\e < 9^^^) is given by 
O with L{e) = Lmax X min[l,(e'/6ln,i„)-^"s'] where L^^x = 
^'inin(^min/ ^max)~^"^' ■ A similar dispersion is introduced in the 
UJ model. The predicted logA^-logS' distribution depends on 
the values of Ausj, (^c, ^'max, f and the normalization O^^^'Lq, 
where the last parameter is determined through a fit to obser- 
vations and is not considered to be a free parameter The scat- 
ter a is constrained by observations, (1.0, 0.9, 0.8, 0.9) J; cr < 
(1.2, 1.2, 1.3, 1.5) for Ausj = (1.5, 2, 2.5, 3), but was allowed 
to vary over a wider range when performing the fit to the data. 
The smallest observed values of 6 provide an upper limit on 
of Oc < 0.05 while the largest observed values of 9 provide 
a lower limit on 6'max of ^max ^ 0.5. For the USJ model the 
shape of the afterglow light curves implies 1.5 < ausj 2-5 
and therefore 1 .5 < Ausj ^ 2.5, however, since this limit does 
not apply to the UJ model and since we wanted to properly 
check the consistency of the USJ with the data we allowed 
Ausj to vary over a wider range when performing the fit to the 
data. 

Combining the constraint from the Frail relation (Eq. |3} 
with the limit on A from the luminosity function allows us 
separately to constrain a and b. From Eq. (|3} we have 
a = 2(4-/fe)-(3-/fe)A = 8-3Afor;t = and 4- A for = 2. 
Since for the USJ model 1 .5 < ausj ^2.5 from the afterglow 
lightcurves, this implies 11/6 < Ausj ^13/6 for k = Q and 
1.5 < Ausj S 2-5 for k = 2. For A; = 0, A must be close to 2 
and b = (4-k)(X-2) falls in the range -2/3<b< 2/3. for 
k = 2 this implies -I ^ b < I. As discussed in |2] direct es- 
timates of e-y suggest < b < I (this applies both to the USJ 
model and to the UJ model). The different constraints on the 
parameters a and b are summarized in Fig. |5] 

If the luminosity function forces Ausj to values less than 2, 
then b must be negative and if Ausj approaches 1.5 this would 
favor an external medium with a k = 2 density profile. Note 
that b < corresponds to the dissipation efficiency increas- 
ing with angle, as would be expected if it is associated with 
a shear layer well outside the jet core. We note that highly 
relativistic jets passing through an external pressure gradient 



less steep than r"'* can develop a shocked layer near the jet 
boundary, due to the loss of causal contact between the jet in- 
terior and wall. If this shock is responsible for the gamma-ray 
emission, it could lead to radiative efficiency increasing with 
9. 

On the other hand, if the luminosity function implies Ausj > 
2 this would imply b > 0, i.e. a gamma-ray efficiency that 
decreases with 9, which is more consistent with direct esti- 
mates of e^. If Ausj approaches 2.5 this would favor a A; = 2 
density profile. 

We performed two fits to the data using a single power-law 
luminosity function which can arise either in the USJ model 
or in the UJ model. In the first fit we let all four free param- 
eters vary: Ausj, Oc, 6'n,ax and a. In the second fit we held the 
value of A fixed at Ausj = 2 and allowed the remaining three 
parameters to vary: 9c, f^max and cr. The second fit was per- 
formed since a value of Ausj = 2 is expected in the simplest 
version of the USJ model where the gamma-ray efficiency e-^ 
is constant (feusj = 0) and therefore ausj = Ausj = 2 because of 
the Frail relation (Eq. |3}. Therefore, it is interesting to test 
whether this simplest version of the the USJ model is consis- 
tent with the observed logA^-log5 distribution. 

In order to assign a value to the fit we divided the 2204 
GRBs in the GUSBAD catalog into 14 bins according to their 
value of S, the peak photon flux. The first 1 1 bins are equally 
spaced in log 5. In the remaining 3 bins, which correspond to 
the highest values of S, we chose a larger range of S values 
so as to have at least A^min = 40 GRBs in each bin, in order to 
have reasonable Poisson statistics. Since the overall normal- 
ization is an additional free parameter in our fits, the number 
of degrees of freedom (d.o.f) in our fits is ly = 14 - (4 + 1) = 9 
in our first fit, and 1^= 10 in our second fit. 

The results of the fits are presented in Table|4]and in Figures 
121 and |3] When A is free to vary we obtain a best fit value of 
Ausj = 2.9;!;o 5, however, when we fix Ausj = 2 we still get an 
acceptable fit. This suggests that although values of Ausj > 2 
are preferred by the data, a value of Ausj = 2 is still possible. 
On the other hand, values of Ausj < 2 become increasingly 
hard to reconcile with the data. Values of Ausj > 2, which 
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0.5 1 1.5 2 2.5 3 

Fig. 2. — The observed cumulative logjQA'-log[Q(Cmax/C[nin) distribution 
taken from the GUSBAD catalog (solid step-like line) compared to the pre- 
dicted log[QA'-log|Q(P/Pii„,) distributions for our best fit models: a single 
power-law luminosity function (relevant for both the USJ and the UJ models), 
both when letting all parameters vary (dashes line) and when fixing Ausj = 2 
(dotted line), and a broken power-law luminosity function (solid line; relevant 
only for the UJ model). 



0.5 1 1.5 2 2.5 

Fig. 3. — The same as Figure |2] but for the differential distribution: 
log|Q(P/Piin,) is divided into 14 bins, the first 11 of equal size, and the re- 
maining 3 with varying sizes chosen such that the number of busts per bin in 
the observed sample is at least W^jn = 40 (in order to have reasonable statis- 
tics so that the Poisson error will not to be too large), and N is the observed 
(circles) or theoretical (lines; only the values at the center of each bin count) 
number of bursts in each bin. The edges of the bins are also plotted. 



are preferred by the data, correspond to > (a gamma-ray 
efficiency that decreases with 6) and for the USJ model 
Ausj ~ 2.5 favors k = 2 (i.e., a stellar wind external density 
profile). 

The first fit gives a = O.S^I'qI which indicates that this fit is 
not very sensitive to the value of a. It also includes the range 
0.8 < fj < 1.3 for Ausj = 2.5 and 0.9 < a < 1.5 for Ausj = 
2.9 that are implied by observations. In contrast, the second 
fit with a fixed Ausj = 2 gives a = 0.8;!;q2 which requires a 
positive value of a in order to get a reasonable fit to the data. 

The best fit values of dc and 6'niax are slightly higher for the 
first fit, but are rather close between the two fits. The best fit 
value of 9c = 0.05-^-^^ for a free A and 0,. = 0.03!{J:[Jf for a 
fixed Ausj = 2 are close to the upper limit of 9c < 0.05 from 
observations, and the allowed confidence regions do not allow 
values much smaller (by more than a factor of ^ 2) than this 
limit. The best fit values of 6'n,ax = 0.7 ± 0.2 for a free A and 
^max = 0.5 ± 0.05 for a fixed Ausj = 2 are close to the lower 
limit of (?inax ^ 0.5 from observations and are not consistent 
with 6lmax = 7r/2. 

The true GRB rate that is implied from our fits is rather 
close to the value of /?grb(z = 0) ^ 0.5 Gpc""* yr"' that was 
found by Perna. Sari & Frail (2003). We obtain a slightly 
larger rate for the free A fit compared to the fixed Ausj = 2 
fit, but the difference is very small. The rates we obtain for 
the USJ m odel are low er by a factor of ^ 300 compared to the 
estimate of Frai l et'an(l2001l) for the UJ model. 

Finally, since the exact choice for the number and sizes 
of the bins used in our fit is somewhat arbitrary and might 
have some effect on our results, we estimate this effect more 
quantitatively by repeating our fit for a single power law lu- 
minosity function with a larger number of bins, 26 instead of 
14 (see Fig. |^. The best fit parameter values remained the 
same as for the original binning (14 bins) up to the significant 
digits shown in Table |4] while the 1 a confidence intervals 
sUghtly changed: Ausj = 2.9!° 2, cr = 0.5^°;^, 9c = 0.05 ± 0.05 
and fi^max = OJto j- The GRB rate is slightly higher, Rgrb{z = 
0) = 0.95;!;|j 22 Gpc"^ yr"', but well within the 1 a confidence 
interval for the original binning. This demonstrates that while 



0.5 1 1.5 2 2.5 

Fig. 4. — The same as Figurel3lbut for a different binning of log[Q(P/Pii„,): 
26 bins and N^i^ =15. 



the exact choice of binning has some effect on our results, this 
effect is rather small. 

In order to estimate the effect of our choice for the star for- 
mation rate (or more accurately the GRB rate that is assumed 
to follow the star formation rate up to an unknown normal- 
ization constant) we repeated our fit with the original bin- 
ning ( 14 bins) bu t with a different star formation rate: SFR2 
from Porciani & Madau ( 2001) . We obtained: Ausj = 2.9![i^, 
cr = 0.5 ± 0.5, 9c = 0.04 ± 0.005 and 6'max = 0.9 ± 0.2. While 
changing the star formation rate has a larger effect than chang- 
ing the binning, the parameter values for the two different star 
formation rates that we considered are still consistent with 
each other within their 1 a confidence intervals. For the new 
star formation rate we obtain a slightly lower true GRB rate: 
Rcrb(z = 0) = 0.65!":?? Gpc-^^ yr-i. 

5. THE LUMINOSITY FUNCTION FOR THE UNIFORM JET MODEL 

The two fits that were discussed in |3 for the USJ model 
still apply to the UJ model with 9c S^in and Ausj 
Auj X 2/(3- q). The second substitution indicates a degen- 
eracy where the luminosity function P{L) and therefore the 
log/V-logS' distribution would be the same for UJ models 
with the same value of Xuj/{3-q) = {auj + bu])/{3-q). This 
degeneracy does not enable us to determine P{9j), which is 
parameterized by q from a fit to the observed logA^-log5 
distribution. Such a fit can only determine rj = 1 + 2/ Ausj = 
1 + (3 - ^) /Auj where P{L) cx L"'' . The degeneracy might be 
broken if we could estimate the true GRB rate in an indepen- 
dent way. This is since the true GRB rate for the USJ model is 
determined from the fit to the observed logA^-logS' distribu- 
tion, and the correction factor for the UJ model depends only 
on q (see Eq. ^^and Figure^. Thus an independent estimate 
of the true GRB rate would enable the determination of q and 
would therefore constrain Pi9j). 

The value of X = a+b is the same for the UJ and USJ models 
if and only if ^ = 1 . In fact, the Frail relation as manifested 
in Eq. ^ implies that for q = I and a given value of k, it 
is enough that either a or b he the same for the two models 
in order for A = a + fe to be the same. Add to this the fact 
that for q = 1 there is an equal number of uniform jets per 
logarithmic interval in 9j, and q = I might be considered as 
the most 'natural' value for q. However, we emphasize that 
there is no physical basis for comparing pairs of models with 
Ausj = Auj, and therefore no a prior reason for choosing q=l. 
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Fig . 5 . — A summary of the constraints on the power law indexes a and b of 
the kinetic energy per solid angle and the gamma-ray efficiency, respectively: 
e oc 9^" and t-y oc fi"*, where 9 = 9gi,s for the US J model and 9 = 9j for the 
UJ model. The constraint from the Frail et al. relation (Eq. |3) is shown 
both for a constant density external medium {k = 0) and for a stellar wind 
environment (k = 2). Also shown are the values of A = a + fc for the best fit 
value (Ausj = 2.9) and the 1 cr confidence interval (2.4 < A^sj < 3.1) for the 
USJ model, as well as A = 2 which is still acceptable. The shaded regions 
show the constraints 1 .5 < ousj ^2.5 from the afterglow light curves, which 
applies only to the USJ model, and < < 1 from estimates of the gamma- 
ray efficiency, which applies to both the USJ model and the UJ model. The 
intersection of these two shaded regions is indicated. 



The only meaningful comparison is to observations. 

The discussion in |4]about the values of Ausj, cr, Omux and 
9c for the USJ model are still valid for the UJ model, where 
9c is replaced by 6'nim and Ausj is replaced by 2Auj/(3-g'). 
The true GRB rate for the UJ model is higher than that for 
the USJ model for the same fit to the observed log A^- log 5 
distribution by a factor of f(q) which is given in Eq. (I14> and 
shown in Figure [2 In Table |4] we provide the values for q = I 
as an example. The true GRB rate that we obtain for ^7 = 1 is a 
factor of « 6-6.6 lower than that of Guetta. Piran & Waxman 
(12004') and a factor of w 45 - 50 smaller than that of Fra il et al.l 
(|^01). However, since the GRB rate that is obtained in this 
way has a strong dependence on the value of q (see Figure 
Q), the rates for q = I are not very meaningful (as discussed 
in the previous paragraph). A constant efficiency (b = 0) and 
a constant energy (a = 2) imply Auj = 2 which, together with 
the best fit value of Ausj ~ 2.9, give q = 3-2Auj/Ausj ~ 1-6. 
For the best fit parameter values this would imply f{q) w 14.4 
and a true GRB rate of /?grb(z = 0) w 12.3 Gpc"^ yr"'. 

It is remarkable that we obtain a good fit to the data for 
the UJ model for a single power-law distribution of jet half- 
opening angles, P{6j), and a scatter a in the peak luminosity 
for a given 9j which is consistent with observations. Previ- 
ous works used a broken power-law luminosity function P{L) 
which corresponds to a broken power-law Pifij), and no scat- 
ter (cr = 0). For this reason, we also performed a fit with 
(T = and a broken power pow er-law lum inosity function P(L) 
with the parameterization of lGuetta. Piran & Waxmanlll2004h . 
where'^ p{L) cx L"'"" for L < L* and P{L) cx L"'"'' for L > L„ . 
The results of this fit are shown in Table0]and Figures|2]and 
We get a good fit where a = 0.6 ± 0.1 and /3 = 0.8!|J;? are 
consistent with having the same value, a « /3 « 0.7. This last 
value implies 77 « 1.7, which is consistent with our best fit 
value of Ausj = 2.9 that implies r/ = 1 +2/Ausj ~ 1.7. This 

*^ Here we define a and fi with a minus sign with respect to their definition 
in lGuetta. Piran & Waxnia <200i) . 



suggests that a broken power-law is not really necessary, and 
a single power-law can provide a reasonable fit to the data 
even for ct = (as we obtain for the first fit from the previ- 
ous section, with different Lmin and Lmax which correspond to 
different ^^ax and 9c). 

6. DISCUSSION 

We have developed a formalism that enables us to cal- 
culate the theoretical logA^-logS' distribution for any GRB 
jet structure, and directly provides the true GRB rate from 
a fit to the observed logA^-logS' distribution. This is a 
more straightforward approach compared to pre vious works 
JGuetta. Piran & Waxmanll2004t lSchnudtll200l which first 
calculated the GRB rate assuming an isotropic emission and 
then introduced a 'correction factor' in order to account for 
the effects of the GRB jet structure. We have applied this 
formalism to the uniform jet (UJ) model and to the universal 
structured jet (USJ) model and performed fits to the GUSBAD 
catalog which includes 2204 BATSE bursts. Our analysis im- 
proves on previous works by: (1) allowing the efficiency in 
producing 7-rays to vary with 9, (2) including an internal dis- 
persion in the peak luminosity L at any given 9, (3) introduc- 
ing an inner core angle 9c and an outer edge 9^^^ in the USJ 
model, and (4) using a larger and more uniform GRB sample. 
The results of our fits are summarized in Table|4]and Figures 
ElandE 

A power law luminosity function has also been fit- 
ted to the the logA^-log5 distribution in some previ- 
ous works jH akkila et al. 1996; Loredo & Wasserman 199^ 
iStern. Tikhomirova & Sven sson 2002, hereafter STS02). Of 
these works, it is most useful to compare our results to those 
of STS02, since their assumptions are the closest to ours: sim- 
ilar to us, they assumed (f^Mii^A) = (0.3,0.7) and that the 
GRB rate follows the star formation rate. It is important to 
keep in mind, however, that STS02 used a different GRB sam- 
ple. STS02 obtained a power law luminosity function at low 
luminosities, with an index of 77 = 1 .4, where P(L) cx L"''. This 
is rather close to the power law index of 77 = 1.7 (or between 
1.6 and 1.8 at 1 a) that we obtain over the whole luminos- 
ity range for a single power law luminosity function. It is also 
close to the power law index of = 1 .6 (or between 1 .5 and 1 .7 
at 1 ct) that we obtain for low luminosities when using a bro- 
ken power law luminosity function. Also, for a single power 
law, we have a maximal luminosity (corresponding to the core 
angle 9c in the USJ model or the minimal jet half-opening an- 
gle 6'n,in in the UJ model) with and exponential tail [due to the 
assumed scatter around the mean value of Lifi)}, which is not 
very different from an exponential cutoff (or decline) at high 
luminosities, that was found to be a viable option by STS02. 
Thus, we conclude that our results are consistent with those 
ofSTS02. 

The main differences between our work and previous works 
in the literature which aimed constraining the GRB luminosity 
function using the peak flux distribution are: (i) the different 
sample that we use, and (ii) the fact that we consider a differ- 
ential peak flux distribution instead of the cumulative distri- 
bution, since random errors propagate in an unknown way in 
the cumulative distribution. Moreover, differently from pre- 
vious authors, we consider the possibility that the jet is struc- 
tured and derive the corresponding luminosity function which 
in turn puts constraints on the jet structure through the fit to 
the logA^-logS' distribution. 

The values we obtain for 9c in the USJ model or 6'inin in the 
UJ model are close to the upper limit of ^ 0.05 that is implied 
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by the smallest values of 6 that are inferred from afterglow 
observations. Furthermore, they are not consistent with zero. 
The value of 6'niax that we obtain is close to the lower limit of 
^ 0.5 from afterglow observations, and is not consistent with 
7r/2. Therefore, we find that including 9c (or 6'niin) and ^?max is 
required in order to obtain a good fit to the data. 

We fit the observed distribution of 9^L(d) (where L is the 
peak luminosity), for a sample of 19 GRBs with a known 
redshift z and an estimate for 6, to a log-normal distribu- 
tion with a standard deviation a. We found the fit to a 
log-normal distribution to be acceptable. The implied val- 
ues of cr are (1.0,0.9,0.8,0.9) < a < (1.2, 1.2, 1.3, 1.5) for 
A = (1 .5, 2, 2.5, 3), where we have taken into account that part 
of the observed dispersion might arise due to error in the es- 
timated value of 6. This is consistent with the wide range of 
possible a values we obtain in our fit for a free A, cr = 0.5;!;q g, 
and with the smaller range of values for our fit with a fixed 
Aim=2,<T = 0.8+ "t 



The Frai fetaD (12001 relation constrains the values of the 
power law indexes a and b of the energy per solid angle in the 
jet (e) and the gamma-ray efficiency (e^), respectively (see Eq. 
|3}- For the USJ model the observed shapes of the afterglow 
light curves imply 1 .5 < ausj < 2.5 jr.ranot & KumaHl20?)3t 
[Kumar & Granot 2003), which in turn implies 1 1/6 < Ausj ^ 
13/6 and -2 /3 < ^usj < 2/3 for ;t = or 1 .5 < Ausj < 2.5 and 
~1 S ^usj S 1 for k = 2. As discussed in ^ direct estimates 
of e.y from afterglow observations suggest < ^ < 1. Our best 
fit value of Ausj = favor a stellar wind external density 

profile (k = 2) and an efficiency that decreases with 9 (b> 0). 
However, Ausj = 2, which corresponds to b = and does not 
constrain the value of k, still provides an acceptable fit to the 
data. The constraints on a and b are summarized in Fig. |5] 

For the UJ model we find a degeneracy in the luminosity 
function: P(L) oc where 77 = l + (3-q)/ Auj and Auj = 
fluj + ^uj- Since a fit to the observed logA^-log5 distribu- 
tion can only constrain the luminosity function P(L), in our 
case it can only provide the value of r/. However, this con- 
strains only the value of (3-^)/Auj = i3-q)/(ai]j + b\jj), and 
therefore still does not enable us to determine the values q and 
Auj separately. An independent estimate of the true ORB rate 
would constrain f{q) and therefore q, and would thus enable 



one to break this degeneracy. 

Alternatively, one can make an additional assumption on 
fluj or bus. The true energy in the jet scales as £ cx 0^~"^' 
and a constant energy corresponds to auj = 2. If we as- 
sume a constant energy, the Frail relation (Eq. |3} implies 
a constant efficiency (bus = 0), and vice versa. In this case 
Auj = 2 and the best fit value of r/ = 1 +2/ Ausj « 1.7 im- 
plies q Ki 1.6. For the best fit parameter values this would 
imply f(q) w 14.4 and a true GRB rate of /?grb(z = 0) w 
12.3 Gpc"^ yr"'. This last value is a factor of « 2.7 lower 
than that of iGuetta. Piran & Waxma nl(l200l and a factor of 
w 20 smaller than that of Frail et al.l 11200 ijT It is also a fac- 
tor of « 4.9 X lO-' smaller than the rate of SNe Type Ib/c, 
^SNib/c(z = 0)«6xl04Gpc-3 yr-1. 

For ouj = 2 and b\js = 0, the 1 a confidence interval of 
2.4 < Ausj < 3.1 in Ausj corresponds to 1.3 < < 1.7, 9.5 < 
/(^)< 16.5,and/;GRB(z = 0)^(7.7-16.5)Gpc-^^yr-i. These 
values of q are not very far from q = I for which there is 
an equal number of jets per logarithmic interval in 9j. It 
implies that there are more jets with Oj ^ ^min compared to 
Oj ~ 6'max, by a factor of (6'max/^'min)^~' 2.1 -9.1 where 

9 < 0rnax/^m,n < 23 (sCC TablcHj. If we allow the value of ei- 
ther fluj or buj to vary, they must both vary together in order to 
satisfy the Frail relation (Eq. |3}- This would cause the values 
of q, f(q) and the true GRB rate for the UJ model to all vary 
accordingly. The strong dependence of f(q) on q (see Figure 
n) implies that this can potentially increase the true GRB rate 
by a large factor. The fact that f(q) > 1 for all q implies that 
the true GRB rate for the USJ model, which is given in Table 
13 provides a lower limit for the UJ model. 



We thank the referee for useful comments which helped im- 
prove the paper. D. G. thanks the Weizmann institute for the 
pleasant hospitality. This work was supported by the RTN 
"GRBs - Enigma and a Tool" and by an ISF grant for a Is- 
raeli Center for High Energy Astrophysics (D. G.), by the 
W. M. Keck foundation and by US Department of Energy 
under contract DE-AC03-76SF00515 (J. G.) and NSF grants 
AST-0307502 (D. G. and M. C. B.) and PHY-0070928 (J. G.). 



REFERENCES 



Bloom, J. S., Frail, D. A., & Kulkami, S. R. 2003, ApJ, 594, 674 

Firmani, C, Avila-Reese, V., Ghisellini, G., & Tutukov, A. V. 2004, submitted 

to ApJ ( astro-ph/0401580 i 
Frail, D.A., et al. 2001, ApJ, 562, L55 

Granot, J., et al. 2001, in GRBs in the Afterglow Era, ed. E. Costa, F. Frontera, 

& J. Hjorth (Berlin: Springer), 312 
Granot, J., Panaitescu, A., Kumar, R, & Woosley, S.E. 2002, ApJ, 570, L61 
Granot, J., & Konigl, A. 2003, ApJ, 594, L83 
Granot, J., & Kumar, R, 2003, ApJ, 591, 1086 

Guetta, D., Piran, T., & Waxman, E. 2004, ApJ, in press lastro-ph/03 11488} 

Guetta, D., Spada, M., & Waxman, E. 2001, ApJ, 557, 395 

Hakkila, J., et al. 1996, ApJ, 462, 125 

Kouveliotou, C, et al. 1993, ApJ, 431, 101 

Kumar, R, & Granot, J. 2003, 591, 1075 

Kumar, R, & Panaitescu, A. 2000, ApJ, 541, L9 

Lipunov, V. M., Postnov, K. A., & Prokhorov, M. E. 2001, Astron. Rep., 45, 

236 

Lloyd-Ronning, N. M., & Zhang, B. 2004, ApJ, in press |astro-ph/0404107 i 
Loredo, T. J., & Wasserman, I. M. 1998, ApJ, 502, 75 
Moderski, R., Sikora, M., & Bulik, T. 2000. ApJ, 529, 151 
Nakar, E., Granot, J., & Guetta, D. 2004, ApJ, 606, L37 
Paciesas, W., et al. 1999, ApJS, 122, 465 



Panaitescu, A., & Kumar, R 2001, ApJ, 554, 667 
Panaitescu, A., & Kumar, R 2002, ApJ, 571, 779 
Panaitescu, A., & Mt'eszf aros, R 1999, ApJ, 526, 707 
Perna, R., Sari, R. & Frail, D., 2003, ApJ, 594, 379 
Porciani, C, & Madau, R 2001, ApJ, 548, 522 
Rhoads, J. E., 1997, ApJ. 487, LI 
Rhoads, J. E., 1999, ApJ, 525. 737 

Rossi, E., Lazzati, D., & Rees, M. J. 2002, MNRAS, 332, 945 

Rossi, E., Lazzati, D., Salmonson, J. D., & Ghisellini, G. 2004, submitted to 

MNRAS I astro-ph/0401 124 1 
Rowan-Robinson, M. 1999, Ap&SS, 266, 291 
Sari, R. 1997, ApJ, 489, L37 

Sari, R., Piran, T. & Halpem, J. 1999, ApJ, 519, L17 
Schmidt, M. 2001, ApJ, 552, 36 

Schmidt, M. 2004, submitted to ApJ I astro-ph/0406^T^ 

Sethi, S., & Bhargavi, S.G. 2001, A&A, 376, 10 

Stern. B. E., Tikhomirova, Ya., & Svensson, R. 2002, ApJ, 573, 75 

Van Putten, M., & Regimbau, T. 2003, ApJ, 593, L15 

Zhang, B., & Meszaros, R 2002, ApJ, 571, 876 



